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Abstract 

i-C | Public opinion is often affected by the presence of committed groups of individuals dedicated to competing 

points of view. Using a model of pairwise social influence, we study how the presence of such groups 
within social networks affects the outcome and the speed of evolution of the overall opinion on the network. 
Earlier work indicated that a single committed group within a dense social network can cause the entire 
network to quickly adopt the group's opinion (in times scaling logarithmically with the network size), so 
long as the committed group constitutes more than about 10% of the population (with the findings being 
qualitatively similar for sparse networks as well) . Here we study the more general case of opinion evolution 
^Lj" when two groups committed to distinct, competing opinions A and B, and constituting fractions pa and 

Pb of the total population respectively, are present in the network. We show for stylized social networks 
£U (including Erdds-Rcnyi random graphs and Barabasi- Albert scale-free networks) that the phase diagram 

of this system in parameter space (paiPb) consists of two regions, one where two stable steady-states 
coexist, and the remaining where only a single stable steady-state exists. These two regions are separated 
by two fold-bifurcation (spinodal) lines which meet tangcntially and terminate at a cusp (critical point). 
We provide further insights to the phase diagram and to the nature of the underlying phase transitions by 
investigating the model on infinite (mean- field limit), finite complete graphs and finite sparse networks. 
For the latter case, we also derive the scaling exponent associated with the exponential growth of switching 
times as a function of the distance from the critical point. 

Introduction 

Since the seminal work of Gabriel Tardc [Tj in the late 1800s, the shaping of public opinion through 
interpersonal influence and conformity has been a subject of significant interest in sociology. This topic 
is especially relevant today due to the preponderance of online social media where individuals can in- 
fluence and be influenced by their numerous and geographically scattered contacts. Public opinion on 
an issue is often shaped by the actions of groups that rigidly advocate competing points of view. The 
most evident example of such a process occurs during elections when multiple parties campaign to influ- 
ence and win over the majority of voters. In this as well as other common scenarios, the predominant 
means of influencing public opinion involves some form of broadcast outreach such as television adver- 
tising, public demonstrations etc. However, even though factors exogenous to the network may have a 
significant effect on individuals becoming informed and engaged in particular issues [2], there is reason 
to believe that large scale changes in behavior or opinion are driven primarily through interpersonal 
influence events occurring within the network. Specifically in the context of rural campaigns, there is 
evidence that interpersonal channels constitute the dominant pathways for effecting individual behavior 
change, even when direct external influence is present [3]. Furthermore, with data on social networks 
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becoming increasingly accessible, there has been a surge of interest in understanding how campaigns can 
be successfully won by leveraging pathways of social influence within the network, thus diminishing the 
need for, or complementing the effect of broadcast outreach. 

Motivated by these observations, we study a simple model that enables us to draw useful insights 
on the evolution of opinions on a social network in the presence of two groups within the network that 
are committed to distinct, competing opinions on an issue. Within the limits of our model, one of the 
questions our work answers is the following. Suppose the majority of individuals on a social network 
subscribe to a particular opinion on a given issue, and additionally some fraction of this majority are 
unshakeable in their commitment to the opinion. Then, what should be the minimal fractional size of 
a competing committed group in order to effect a fast reversal in the majority opinion? In addition to 
answering this question quantitatively, we show the existence of two distinct types of phase transitions 
that can occur in the space of committed fraction pair values. 

We model the dynamics of social influence using a two-opinion variant of the Naming Game [H [5] [5] 
which also corresponds to a special case of the game introduced and studied in [TJ [5] . The same model 
was referred to as the binary-agreement model in [9]. In this model, at any time, a node possesses cither 
one of the two competing opinions (i.e. the node is in state A or state B), or both opinions simultaneously 
(state AB). In a given time step, we choose a node randomly, designate it as the speaker and choose one 
of its neighbors randomly and designate it as the listener. The speaker proceeds to convey its opinion 
to the listener (chosen randomly if it possesses two) to the listener. If the listener possesses this opinion 
already, both speaker and listener retain it while eliminating all other opinions; otherwise, the listener 
adds the opinion to his list. A table of possible interactions and outcomes between node-pairs is provided 
in Table SI. We emphasize that each node interacts and is influenced only by its neighbors on the 
network. There is no element in our model that represents an external influence mechanism such as the 
use of media, public demonstrations, or door-to-door campaigns by members of the competing groups. 
Except for their being un-influencable, the committed nodes are assumed to be identical in all other 
respects to uncommitted nodes. In particular, committed nodes do not influence their neighbors at a 
different rate or with a higher strength than uncommitted nodes. 

Opinion dynamics models involving committed individuals all subscribing to a unique opinion have 
been studied previously in [9j [lOj [11] [12] . The situation pertinent to this paper - that of two competing 
committed groups - has received considerably greater attention (TTJ [T21 HH HH US] . Mobilia et al. [33] 
studied how the presence of zealots (equivalent to committed individuals) affected the eventual distribu- 
tion of opinions (stationary magnetization) in the case of the voter model. They demonstrated that the 
distribution for a finite sized network was Gaussian, with a width inversely proportional to the square 
root of the number of zealots, and centered at * + T*~ where z + ,z_, represent the fraction of zealots in the 
two competing states. Similarly to [H], Yildiz et al. [16] studied the properties of steady-state opinion 
distribution for the voter model with stubborn agents, but additionally considered the optimal placement 
of stubborn agents so as to maximally affect the steady-state opinion on the network. Interestingly, unlike 
in the model studied here, in the voter model, no transitions in steady-state magnetization are observed 
as the committed fraction pair values are smoothly varied. Biswas et al. [15] considered the effect of 
having rigid individuals in a one-dimensional system of binary opinion evolution, and demonstrated a 
power-law dependence for the decay of steady-state magnetization on the fraction of rigid individuals. The 
work done in jll| 113] is similar in spirit to our work here; however, an important difference is that these 
studies only considered the infinite-network size limit for complete graphs. We study finite networks, both 
complete and sparse, and provide semi-analytical arguments regarding timcscalcs that become relevant 
when the network size is finite. 
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Analysis 

First, we study the mean-field version of the model, also being equivalent to the dynamics on the complete 
graph in the limit of infinite system size. We designate the densities of uncommitted agents in the states 
A, B and AB by ua, n B and uab- We also designate the fraction of nodes committed to state A, B 
by pa, Pb respectively. These quantities naturally obey the condition: ua + n B + tlab + Pa + Pb = 1- 
In the asymptotic limit of network size, and neglecting fluctuations and correlations, the system can be 
described by the following mean- field equations, for given values of the parameters pa and p B : 

2 3 
n A n B + n AB + n A n AB + A n ab ~ PbUa 

2 3 
n A n B + n AB + n B n A B + -PbTIab - PA n B 

(1) 

The evolution of uab follows from the constraint on densities defined above. In general, the evolution of 
the system depends on the relative values of pa and p B . In the case of pa > 0, ps = (or equivalcntly, 
Pb > 0, pa = 0) there is only a single group of committed nodes in the network, all of whom subscribe to 
the same opinion. This was the case studied in [IT] [9] [12]. In this scenario, a transition is observed when 
this committed group constitutes a critical fraction of the total network. Specifically, the transition point 
separates two dynamical scenarios in the phase space, (tiajTIb), of uncommitted node densities. Below 
the critical value, the absorbing state (e.g., ua = 1 — Pa, n B = uab = when pa > 0,Pb = 0) coexists 
in phase space with a stable mixed steady-state and an unstable fixed ("saddle") point. At or above the 
critical value, the latter non-absorbing steady-state and the saddle point cease to exist. Consequently, 
for a finite system, reaching the (all A) consensus state requires an exponentially long time when p is 
less than the critical value. Beyond the critical value this time grows only logarithmically with network 
size. Note that this critical value or threshold is analogous to a spinodal point [XT] [18] associated with 
an underlying first-order (or discontinuous) transition in equilibrium systems. 

In order to effectively characterize the behavior of the system governed by Eqs. (UJ for pa,Pb > 0, we 
systematically explore the parameter space (j>a,Pb) by dividing it into a grid with a resolution of 0.000125 
along each dimension. We then numerically integrate Eqs. (JTJ) for each (pa,Pb) pair on this grid, assuming 
two distinct initial conditions, Ua = 1 — Pa — Pb, ^b = it-ab = and ns = 1 — Pa — PB,nA = tiab = 0, 
representing diagonally opposite extremes in phase space. The results of this procedure reveal the picture 
shown in Fig. [JJ in different regions of parameter space. As is obvious, with non-zero values for both 
Pa,Pb, consensus on a single opinion can never be reached, and therefore all fixed points (steady-states) 
are non-absorbing. With (j>a,Pb) values within the region denoted as I which we refer to as the "beak" 
(borrowing terminology used in [J), the phase space contains two stable fixed points, separated by a 
saddle point, while outside the beak, in region II, only a single stable fixed point exists in phase space. 
In region I, one fixed point corresponds to a state where opinion A is the majority opinion (A-dominant) 
while the other fixed point corresponds to a state where opinion B constitutes the majority opinion 
(5-dominant). Figure [JJ shows representative trajectories and fixed points in phase space, in different 
regions of parameter space. Similar phase diagrams have been found in other two-parameter systems in 
different contexts including chemical reactions [4] and genetic switches [20] . 

In order to study the nature of the transitions that occur when we cross the boundaries of the beak, we 
parametrize the system by denoting ps = cpa where c is a real number. Then, we systematically analyze 
the transitions occurring in two cases: (i) c = 1 and (ii) c ^ 1. It can be shown that along the diagonal 
line c = 1 the system undergoes a cusp bifurcation at pa = Pb = 0.1623. The movement of the fixed 
points as pa and ps are smoothly varied along the diagonal line is shown in Figure SI. Henceforth, we 
denote the value of pa and pb at the cusp as p c . As is well known, at the cusp bifurcation two branches of 
a saddle- node (or fold) bifurcation meet tangentially [TJ . These two bifurcation curves form the boundary 
of the beak shown in Fig. [1] A detailed analysis demonstrating that pa — Pb — Pc constitutes a cusp 
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Figure 1. Mean-field picture in parameter space. The phase diagram obtained by integrating the 
mean- field Eqs. (Q]). The two lines indicate saddle- node bifurcation lines which form the boundary 
between two regions with markedly different behavior in phase space. For any values of parameters 
within the beak, denoted as region I, the system has two stable fixed points separated by a saddle point. 
Outside of the beak, in region II, the system has a single stable fixed point. The saddle-node bifurcation 
lines meet tangentially and terminate at a cusp bifurcation point. 



bifurcation, as well as a semi- analytical derivation of the bifurcation curves is provided in the Supporting 
Text SI (Sections: 1, 2, 3). The cusp bifurcation point is analogous to a second-order (or continuous) 
critical point seen in equilibrium systems, while bifurcation curves are analogous to spinodal transition 
lines. 

Next, we study the stochastic evolution of opinions on finite-sized complete graphs through simula- 
tions. Here, we systematically vary c from 1 to to obtain the right bifurcation curve, and therefore by 
virtue of the A-B symmetry in the system, also obtain the left bifurcation curve. In particular for a given 
value of c we obtain the transition point by varying p A (with ps — cPa) ancl measuring the quantity: 



m = (n B - n A )l {l-p A - p B ) 



(2) 



which we utilize as an order parameter. The above order parameter is analogous to the "magnetization" 
in a spin system as it captures the degree of dominance of opinion B over opinion A and is conventionally 
used to characterize the nature of phase transitions exhibited by such a system. 
Another quantity, the Binder cumulant, defined as 



U N = l 
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Figure 2. Behavior of typical order parameters as a function of linear trajectories of slope 
c that pass through the origin, in parameter space for a complete graph, (a)-(b) Steady-state 
magnetization m denned in the text, for successive pa, Pb pairs along lines of slope c = 1 and c = 0.5 
respectively that pass through the origin. The c = 1 line in parameter space passes through the cusp 
point and gives rise to a second-order phase transition, while the c = 0.5 line passes through a point on 
the (right) bifurcation line giving rise to a first-order phase transition. Here 10 realizations of social 
influence dynamics were performed for each pa, Pb pair, starting from the initial condition 
riA = 0, 71b = 1 — Pa — Pb, and the magnetization was measured conditioned on the system remaining 
in the steady state that it initially converged to. (c)-(d) Scaled variance, Xn, defined in the text for 
successive pa,Pb pairs along lines of slope c = 1 and c = 0.5 respectively, that pass through the origin. 
(e)-(f) Binder cumulant Un defined in the text for successive pa,Pb pairs along lines of slope c = 1 and 
c = 0.5 respectively, that pass through the origin. Data for (c),(d),(e) and (f) were generated from 10 
realizations of the social influence dynamics, per pa, Pb pair, for each of two initial conditions: 
n A = 1 - pa - Pb, n B = and n A — 0, n B = 1 - Pa - Pb- 
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Figure 3. Picture in parameter space for a complete graph obtained from analytical and 
simulation results. The bifurcation lines and the cusp point in parameter space were obtained 
analytically from the mean field equations and are compared with those found using simulations for 
finite-sized complete graphs. Analytical and simulation curves show excellent agreement as N increases. 
The location of the transition occurring across the bifurcation curve was obtained using the Binder 
cumulant Un (Fig- 2(f)), while the location of the cusp point was obtained by using variance of m (Fig. 
2(c)). For both analytical and simulation results, the bifurcation curves arc obtained by identifying the 
critical points that lie on linear trajectories in parameter size described by ps = cpA- This process is 
carried out for different values of c between and 1 at intervals of 0.1, and for each value of c, pa is 
varied at a resolution of 0.002. In simulations, for each such combination of {pa,Pb) obtained, we 
perform averages over 10 realizations of the social influence dynamics, for each of two initial conditions, 
riA = 1 — PA — Pb, and ua = 0, with ub = 1 — fiA — PA — Pb for each case. 

for a system of size N, is commonly used to distinguish between different types of phase transitions |17j . 
The utility of the Binder cumulant comes from the markedly different signatures we expect it to produce 
along a spinodal trajectory (e.g. c = 0.5) - one that passes through the spinodal line - and one along a 
trajectory that passes through the critical point (e.g., along the diagonal, c = 1). This difference arises 
from the following distinction in the evolution of the distribution of m, P(m), along these trajectories. 
Along a spinodal trajectory starting from a point where pa = Pb, an initially symmetric (about m = 0), 
bimodal P(m) becomes asymmetric and unimodal upon crossing the spinodal line, with the single mode 
eventually becoming a delta function. In contrast, along the diagonal trajectory in parameter space, P(m) 
is initially a double-delta distribution (for pa = Pb <C p c ), symmetric about m = 0, and it smoothly 
transitions to a zero-centred gaussian distribution as the critical point is crossed. The definition of Un 
indicates that Un = 2/3 for a delta function distribution (also for a symmetric, double-delta distribution 
about m = 0), while Un — for a zero-centered Gaussian distribution, and thus readily yields the 
limiting Un values at both extremes of the spinodal and diagonal trajectory. As illustrated in Fig. [2j 
Un as a function of pa shows distinct behaviors for c = 1 and c = 0.5, indicating the existence of 
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a second-order (or continuous) transition point at pa = Pb = Pc(N) (Fig. [U(e)) and first-order (or 
discontinuous) phase transition points (Fig. HJf)) along off-diagonal trajectories [T7], respectively. The 
second-order critical point p c (N) converges to the mean-field value, p c sa 0.1623, as N becomes larger. 
The dip observed in Un along the off-diagonal trajectory serves as an excellent estimator of the location 
of the first-order (spinodal) transition for a finite network. Thus, to reiterate, for a finite network, the 
second-order transition point and the first-order transition (spinodal) lines are respective analogues of 
the cusp bifurcation point and the saddle-node bifurcation curves observed in the mean-field case. 

The fluctuations of the quantity m can also be used to identify a transition point, particularly for the 
case of the second-order transition. In particular, in formal analogy with methods employed in the study 
of equilibrium spin systems, the scaled variance: 

X N = N{(\m\-{\m\)) 2 ) (4) 

serves as an excellent estimate for the second-order transition point p c for a finite network. As shown 
in Fig. [He), Xn peaks at a particular value of pa, with the size of the peak growing with N (and 
expected to diverge as N — >• oo). In the case of the spinodal transition, one studies fluctuations of m 
(Xjv = N((m — (m)) 2 )) restricted to the metastable state [22j [23] until the spinodal point (Fig. EJd)) 
at which the metastable state disappears, and fluctuations of m in the unique stable state beyond the 
spinodal point (Fig. [2jd)). 

Figure [3] shows the bifurcation (spinodal) lines obtained via simulations of finite complete graphs 
by using the Binder cumulant (Fig. [2ff))to identify the location of the spinodal phase transition, and 
demonstrates that its agreement with the mean-field curves improves as N grows. The cusp points shown 
here are identified in simulations as the locations where Xn reaches its peak value (Fig. Etc)). 

In the region within the beak, the switching time between the co-existing steady-states represents the 
longest time-scale of relevance in the system. The switching time is defined as the time the system takes 
to escape to a distinct co-existing steady-state, after having been trapped in one of the steady-states 
(see Fig. 0|a)). In stochastic systems exhibiting multistability or metastability, it is well known that 
switching times increase exponentially with N for large N (the weak-noise limit) [6][25j[26l[4] Furthermore, 
the exponential growth rate of the switching time in such cases can be determined using the eikonal 
approximation [4j[27]. The basic idea in the approximation involves (i) assuming an eikonal form for the 
probability of occupying a state far from the steady-state and (ii) smoothness of transition probabilities 
in the master equation of the system. This allows the interpretation of fluctuational trajectories as 
paths conforming to an auxilliary Hamilton- Jacobi system. This in turn enables us to calculate the 
probability of escape allowing an optimal fluctuational path that takes the system from the vicinity of the 
steady-state to the vicinity of the saddle point of the deterministic system. The switching time is simply 
the inverse of the probability of escape along this optimal fluctuational path. We defer details of this 
procedure to Supporting text SI: Section 4. Using this approach we find that for the symmetric case, 
PA = Pb = P < Pc, the exponential growth rate of the switching time s ~ (p c — p) u with v « 1.3 (Fig. [4] 
(c) ). Thus, along the portion of the diagonal within the beak: 

^switching ~ CXp[(p c -p) v N] (5) 

Outside the beak, the time to get arbitrarily close to the sole steady-state value grows logarithmically 
with N (not shown). 

The results presented so far show that there exists a transition in the time needed by a committed 
minority to influence the entire population to adopt its opinion, even in the presence of a committed 
opposition (i.e. in the case where both pa,Pb > 0), as long as pa,Pb < Pc- (Note that the case 
Pa > 0,pb — was considered in [9]). For example, assume that initially all the uncommitted nodes 
adopt opinion B, and that pa = Pb < Pc- Then, the steady-state that the system reaches in ln(iV) time 
is the one in which the majority of nodes hold opinion B. Despite the fact that there exist committed 
agents in state A continuously proselytizing their state, it takes an exponentially long time before a 
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large (spontaneous) fluctuation switches the system to the A-dominant steady-state. For identical initial 
conditions, the picture is qualitatively the same if we increase pa keeping ps fixed, as long as {pa,Pb) 
lies within the beak. However, when (pa,Pb) lies on the bifurcation curve or beyond, the i?-dominant 
steady-state vanishes, and with the same initial conditions - where B is the initial majority - it takes 
the system only In (AT) time to reach the A-dominant state (the only existing steady-state). Thus, for 
every value of an existing committed fraction ps (< p c ) of B nodes, there exists a corresponding critical 
fraction of A nodes beyond which it is guaranteed that the system will reach an A dominant state in 
ln(iV) time, irrespective of the initial conditions. However, for any trajectory in the parameter space 
in a region where either pa or pb is (or both are) greater than p c , no abrupt changes in dominance or 
consensus times are observed. Instead, the dominance of A or B at the single fixed point smoothly varies 
as the associated committed fractions are varied. Moreover, the system always reaches this single fixed 
point in ln(7V) time. 

Finally, we study how opinions evolve in the presence of committed groups on sparse graphs, most 
relevant to social networks. We study Erdos-Renyi (ER) random graphs [25J as well as Barabasi-Albcrt 
networks |29j . For each of these sparse networks, we find the same qualitative behavior as found for 
the complete graph. As shown in Figs. [5j [U as the average degree of the sparse networks increases, the 
bifurcation lines in parameter space tend to approach their mean-field counterparts. Although we do 
not study sparse networks analytically here, we note that in another instance of a phase transition for 
a similar model studied in [TJ, it was demonstrated using heterogeneous mean-field equations that the 
behavior of sparse networks is qualitatively similar to that of complete graphs. Figure [7J visually depicts 
typical instances of the evolution of opinions on an ER random graph for (pa,Pb) values within and 
outside the beak. 



Discussion 

Using a simple model, we have explored and quantified possible outcomes for the evolution of opinions 
on a social network in the presence of groups committed to competing opinions. Broadly speaking, our 
results indicate that as long as the fraction, ps, of nodes committed to a given opinion B is held fixed at a 
value less than a critical value p Cl it is possible to induce the network to quickly tip over to a state where 
it widely adopts a competing opinion A, by introducing a fraction of nodes committed to opinion A. The 
value of the competing committed fraction, pa, at which this tipping point arises depends on the value 
of pb, and is determined by the bifurcation curve (see Fig.[l}. Importantly, for a given value of ps < Pc, 
the excess commitment pa — Pb required for the network to tip over to A is a decreasing function of ps 
that reaches zero when ps = p c - While the critical value p c itself may depend on the network structure 
and its size, the feature described above holds for the three different classes of networks studied here. 
A corollary to this feature is that if the committed fraction ps is held fixed at a value greater than pc, 
increasing the competing committed fraction pa only yields continuous incremental gains in the adoption 
of A (i.e., no tipping point or discontinuous changes in opinions exist). We analytically determine that 
p c = 0.1623 for infinite-sized complete graphs, which as observed from our simulation results in Figs. [5j|6] 
appears to constitute a good upper bound to the value of p c for sparse networks. 

Our results could be of utility in situations where public opinion is deadlocked due to the influence of 
competing committed groups. Perhaps one example of such a situation is the observed lack of consensus 
in the U.S. on the existence of human-induced climate change. Indeed, there is evidence in this particular 
case that the commitment of individuals to particular political ideologies may have an effect on their 
opinions }30j . 

Another scenario to which our model could bear some relevance is the adoption of competing industrial 
standards. Particularly in situations where a network of entities collaborate or are interdependent, there 
is a natural attempt at agreement in standards or protocols between interacting members. A classic 
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Figure 4. Evolution of order parameter m and the exponential growth in switching time 
as a function of distance from the second-order critical point, (a) Switches in the value of m as 
a function of time t for a sample evolution (with initial transient removed) of the system when 
PA = Pb = 0.154 (< p c ). This reflects the system repeatedly switching between the ^-dominant 
steady-state (m > 0) and the B-dominant steady-state (m < 0). (b) Sample evolution of the system 
(with initial transient removed) for pa = Pb = 0.2 (> p c ). The system fluctuates randomly about the 
only existing steady-state in which densities of A and B nodes are equal, (c) The dependence of s in the 
exponential scaling T sw i tc hing ~ exp(sTV) when pa = Pb = P {p < Pc) as a function of (p c — p), obtained 
using the eikonal approximation (see SI: Section 4) 



example of this scenario is the case of the Sellers' screw manufacturing standard that proliferated despite 
competition from the Whitworth standard |31j . A key factor responsible for the eventual success of the 
Sellers standard was William Sellers' leveraging of his connections to corporations and manufacturers |32j . 
whom he persuaded to become adopters of his standard . Furthermore, the network of interdepcndcncics 
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Figure 5. Results for Erdos-Renyi random graphs. (a)Thc bifurcation lines and cusp point in 
parameter space obtained through simulations of Erdos-Renyi random graphs of size N = 5000 with 
different average degrees. The mean-field analytical curve is shown for comparison. For simulation 
results, the bifurcation curves are obtained by identifying the critical points that lie on linear 
trajectories described by ps = cpa in parameter space. This process is carried out for different values of 
c between and 1 at intervals of 0.1, and for each value of c, pa is varied at a resolution of 0.002. For 
each such combination of {jpa,Pb) obtained, we perform averages for quantities of interest over 10 
realizations of networks (with a single realization of the social influence dynamics per network) , for each 
of two initial conditions, ua = 1 — Pa — Pb and ua = with tib = 1 — Ua — Pa — Pb in each case. 
(b)-(c) Steady-state magnetization for ER graphs with (k) = 6 and different sizes N , as parameter pair 
values are varied successively along slope c = 1 and slope c = 0.5 lines in parameter space respectively. 



between industries at the dawn of the mass-manufacturing era played an important role in the adoption of 
the standard becoming widespread. It should be pointed out that in this case, the uncommitted members 
of the population initially adhered to neither standard - this situation can however be accommodated in 
our model by assigning each uncommitted his own unique "opinion" to begin with in close analogy to 
initial conditions for the original Naming Game [4j [5] . 

A more recent example of such a scenario is the competition between Flash and HTML5 in web- 
dcvelopmcnt. There is speculation that Flash, which until recently was the predominant platform for 
animated web content, is gradually ceding its dominance to HTML5 as a result of the increasing market- 
share of Apple's mobile devices which exclusively support the latter |33j . 

A potential competition between DC fast charging standards is also expected as electric vehicles 
become increasingly popular with consumers. The front-runners in the mass manufacture of electric 
vehicles have opted for the CHAdeMO standard, and charging stations compatible with the standard have 
begun proliferating in the US, Europe and Japan [34l [35] . An alternative to CHAdeMO currently being 
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Figure 6. Results for Barabasi- Albert networks. (a)The bifurcation lines and cusp point in 
parameter space obtained through simulations of Barabasi- Albert networks of size TV = 5000 with 
different average degrees. For simulation results, the bifurcation curves are obtained by a similar 
method as described in the legend of Fig. EJa). (b)-(c) Steady-state magnetization for BA networks 
with (k) = 6 and different sizes N , as parameter pair values are varied successively along slope c = 1 
and slope c = 0.5 lines in parameter space respectively. 

developed by the Society for Automotive Engineers (SAE), which governs the development of standards 
in the US automotive industry, is being touted by some car manufacturers as more cost effective as well as 
technologically superior. However, by the time the first cars employing the SAE standard hit the market, 
CHAdeMO charging stations are expected to be rather widespread, thus making a competition between 
the two inevitable (36] [34] . As new collaborations are forged between car-makers especially in the area 
of electric vehicle development [371 [35] ( m addition to collaborations with energy suppliers), the outcome 
of this competition could be significantly influenced by manufacturers who arc already committed to one 
of the two standards through their investment in them. 

To conclude, we have presented results from a simple, abstract model for understanding how opinions 
on a social network evolve through social influence when there are multiple groups within the network 
dedicated to competing opinions. Despite the simplicity of our model, we believe the insights provided 
here form a useful theoretical complement to data-driven studies [39] and randomized evaluations [40] 
aimed at understanding the spread of opinions. 

Note: Subsequent to this paper's initial posting on arxiv.org on 12/31/2011 and its acceptance for 
publication with minor revisions on 1/26/2012, S. Jolad sent us independent unpublished results by 
D. Linford, P. Hochendoner, A. Reagan, and S. Jolad addressing competing committed groups on the 
complete graph. 
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(a) 



(b) 




Figure 7. Visualization of opinion evolutions. The evolution of opinions on an ER random graph 
with N = 200 and (fc) = 6 for two (pa,Pb) pairs. In each case tib = 1 — Pa — Pb and ua — 0. Nodes 
holding opinion A are depicted in red, while nodes holding opinion B arc shown in green. Nodes with 
larger diameters are committed nodes. Top: The case pa = Pb = 0.1 for which the system is in region I 
in parameter space (following the terminology of Fig. [TJ and the system is trapped in a £>-dominant 
steady-state. Even after 2500 time steps, the system continues to remain trapped in this state (inset) 
with riA ~ 0.05. Bottom: The case pa = 0.125 and pb = 0.1 for which the system is in the region II, 
and undergoes an abrupt transition (inset) to the A-dominant state within 250 time steps. 
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1 Analysis of steady states for p^ = p B : existence of a critical 
value p A = p B = p c . 

For notational simplicity we replace ua by x and Ub by y. The mean field equations describing the 
system with pa = Pb = P, < p < 0.5 then are: 

dx 3 

— = -xy + (1 - x - y - 2p) 2 + x(l - x - y - 2p) + -p(l - x - y - 2p) - px 
dx 3 

— = -xy + (l-x-y- 2pf + y(l - x - y - 2p) + -p(l -x-y-2p)-py 

(1) 

where uab = 1 — x — y — 2p. In the steady state, dx/dt = dy/dt = 0, and the resulting equations can 
be solved to yield four solutions for (x,y). Out of these one solution lies outside the valid range for all 
feasible values of p, i.e., < p < 0.5. The valid fixed points for Eqs. [T]are: 



3 i 

2/1 = g _ 2^ 5 ~ 2p ~ P 



1 1 



%2 = ~ 

2 



V2= 2 



2^ 



6p 



1 1 



5^ 



6p 



2 P 



2 P 



x 3 



V3 



V 1 - p 2 - 6p - 
V 1 - p 2 - 6p - 



Since the solutions are symmetric in x and y, in order to investigate the range of p over which these 
solutions are valid, we restrict our analysis to y. The solution y\ is valid for all values of p. For y2, 2/3 to 
be valid solutions, we require U(p) = 1 — p 2 — 6p > 0. U(p) is a monotonically decreasing function for 
p > 0, and the value of p at which U(p) first crosses zero is the critical point. 



Pc 



10- 3 ~ 0.1623. 



(2) 



Thus, there exist three fixed points in the range [0,p c ]. In the range (p c ,0.5] only one valid fixed point 
exists, viz. (x\,yi). 

We can further examine the stability of the obtained fixed points. Linear stability analysis yields the 
following stability matrix 

-2 + 2y 



-1 - £ 
-2 + 2x* + %p 



IP 



(3) 



where is the fixed point under consideration. 
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Before interaction 


After interaction 


A 4 A 


A- A 


A 4 B 


A - AB 


A A AB 


A- A 


B 4 A 


B - AB 


B 4 B 


B - B 


B 4 AB 


B - B 


AB 4 A 


A- A 


AB 4 B 


AB - AB 


AB 4 AB 


A- A 


AB 4 A 


AB - AB 


AB 4 B 


B - B 


AB 4 AB 


B - B 



Supplementary Table 1. Shown here are the possible interactions in the binary agreement model. 
Nodes can possess opinion A, B or AB, and opinion updates occur through repeated selection of 
speaker-listener pairs. Shown in the left column are the opinions of the speaker (first) and listener 
(second) before the interaction, and the opinion voiced by the speaker during the interaction is shown 
above the arrow. The column on right shows the states of the speaker-listener pair after the interaction. 

The eigenvalues of the stability matrix at the fixed point are given by : 

A = -(2+p) ± y / 26p 2 + (20(x* + if) - 36) + 16(1 - x* - y* + x*y*), 

and examination of the real part of these eigenvalues indicates that {x2, 2/2) and (2:3, j/3) are stable fixed 
points, and {x\,y\) is an unstable fixed point (saddle point) for p < p c = 0.1623. For p > p c , (xi,yi), 
the only valid fixed point, is a stable fixed point. Supporting Figure 1 shows the movement of the fixed 
points in the phase space as a function of p. 

2 Existence of a cusp point 

Suppose that a cme-dimcnsional parameter (a) dependent system 

^ = /(x;a),xe3? 1 ,ae5R m (4) 

with smooth function /, has an equilibrium at x — for a = 0, and let f x (0;0) = and f xx (0;0) = 
hold. Further, assume that the non-degeneracy conditions (e.g., f xxx (Q;0) 7^ 0) are satisfied. Then the 
system undergoes a cusp bifurcation at x = [T]. 

We prove that such a cusp bifurcation is encountered in our system (i.e., Eq. [1} at pa = Pb = Pc as 
we move along the diagonal in parameter space {pa = Pb)- Note that our system is two-dimensional. To 
be able to use the above theory, we first need to reduce the dimensionality of our system. The Center 
Manifold Theorem [5] guarantees the existence of a one-dimensional center manifold to which we can 
restrict our system, and such a system preserves the same behavior as the original system in the vicinity 
of the steady-state under consideration. Once we get the restricted system, we can perform the usual 
bifurcation analysis in one-dimensional system. Following this idea, we first shift the coordinates such 
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Supplementary Figure 1. Movement of fixed points as pa and ps are smoothly varied 
along the diagonal line pa = Pb- For pa — Pb,Pc ~ 0.1623 three fixed points exist, two of which are 
stable, and the third is unstable. For pa = Pb > Pc, only a single stable fixed point exists. 



that the origin is located at the critical point we found from the pb = PA case (for simplicity, we denote 
Pa by p and pb by r), i.e., (xq, yo',Po, r o) — (0.2565, 0.2565; \/W— 3, y/lO — 3). In the shifted coordinates, 
the eigenvalues and eigenvectors are given by A = [0; -2.1623] and V = [-0.7071, 0.7071; 0.7071, 0.7071]. 
Using the transformation [x y] T = V[x y] T , and after some algebraic manipulations, we obtain in the 
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dy 
dt 



new co-ordinate system: 
dx 

— = 0.7071(1.5p+0.2434)(p + r+ lAUy- 0.1623) 
-0.7071(1. 5r + 0.2434)(p + r + lAUy - 0.1623) 
-0.7071(0.70715 + 0.7071y + 0.2566)(p + r + lAUy - 0.1623) 
-0.7071(p + 0.1623)(0.70715 + 0.7071y + 0.2566) 
+0.7071(0.7071y - 0.70715 + 0.2566)(p + r + lAUy - 0.1623) 
+0.7071(r + 0.1623)(0.7071y - 0.70715 + 0.2566) 

-0.7071(1. 5r + 0.2434)(p + r + lAUy - 0.1623) 

-0. 7071(1. hp + 0.2434)(p + r + lAUy - 0.1623) 
-0.7071(0.70715 + 0.7071y + 0.2566)(p + r + lAUy - 0.1623) 
-0.7071(p + 0.1623)(0.70715 + 0.7071y + 0.2566) 
-0.7071(0.7071y - 0.70715 + 0.2566)(p + r + lAUy - 0.1623) 
-1.414(0.70715 + 0.7071y + 0.2566)(0.7071y - 0.70715 + 0.2566) 
-0.7071(r + 0.1623)(0.7071y - 0.70715 + 0.2566) 

+ 1.414(p + r + 1.414y- 0.1623) 2 (5) 

Next, we use a quadratic approximation for the center manifold of the above system [2] i.e. we assume 
y = h(x) = \wx 2 . We can find w by comparing two expressions obtained for the first is obtained 
by using ^| = ^| an d then using the first equation in Eq. [5] and the quadratic approximation for y; 
the second is obtained by direct substitution of the quadratic approximation into the second equation in 
Eq. \5\ Doing this yields: 

y = h(x) = -0.70715 2 /(4p + 4r - 2.1620) 
Hence we obtain the following one dimensional system restricted to the one-dimensional center manifold: 



— = 0.1814r-0.1814p-5(1.5p+1.5r) 
dt 

+0.7071(1. 5p + 0.2434)(p + r - 0.1623) 
-0. 7071(1. 5r + 0.2434)(p + r - 0.1623) 

-5 2 (0. 7071(1. 5p + 0.2434)/(4p + 4r - 2.1623) (6) 
-0.7071(1. 5r + 0.2434)/(4p + 4r - 2.1623) 
-0.7071(p+ 0.1623)/(8p + 8r - 4.3246) 
+0.7071(r + 0.1623)/(8p + 8r - 4.3246)) 
+5 3 /(4p + 4r- 2.1623) 

It is easy to check that the origin in this transformed system satisfies the necessary conditions for 
a cusp bifurcation. The origin of this transformed system corresponds to the point pa = Pb = Pc in 
our original system Eq. [T] Thus, the system undergoes a cusp bifurcation at pa — Pb = Pc where 
p c = v/10-3 « 0.1623. 
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3 Mapping out the bifurcation curves (first order transition 
lines) 

In order to map out the first-order transition line (bifurcation curve) we adopt a semi-analytical approach. 
We assume pb = cpa with c < 1 to obtain the lower bifurcation curve (symmetry of the system allows 
us to obtain the upper bifurcation curve, given the lower one). Using Eqs. [TJ the fixed point condition 
becomes (for simplicity, we denote pa by p) : 

f(x,y,p) = -xy + (1 - x - y - (1 + c)p) 2 + x(l - x - y - (1 + c)p) 
3 

+ ~ x — 2/ — (1 + c )p) — C P X = 

g{z,V,p) = ~xy + (1 - x - y - (1 + c)p) 2 + y(l - x - y - (1 + c)p) 
3 

+ o c p(! - x - y - (i + c )p) - py = o 



In addition, for a fold bifurcation, we also require that the stability matrix has an eigenvalue with 
zero real part. Since, the valid solutions in our case are always real, this is equivalent to requiring the 
determinant of the stability matrix to be zero. Thus the condition \Q\ = (with Q given by Eq. [3]) along 
with Eqs. [7| enable us to determine for a given c, the location (pa,cpa) at which the bifurcation occurs. 
By numerically solving these equations for different values of c, < c < 1 at intervals of 0.1, we obtain 
the lower bifurcation curve shown in Fig. 2 of the main text. 



4 Optimal fluctuational paths, the eikonal approximation and 
switching times between co-existing stable states 

The master equation for our system takes the general form: 
dP(X, t) 



at 



E 



W(X - r, r)P(X - r, t) - W(X, r)P(X, t) 



where X = [Na Nb] t denotes the (macro) state of the system as vector whose elements are the numbers 
of uncommitted nodes in state A and B respectively, W(X, r) is the probability of the transition from 
X to X + r, and r runs over the allowed set of displacement vectors in the space of macro-states. For 
our system, r runs over [1 0] T , [0 1] T , [2 0] T , [0 2] T , [—1 0] T , [0 — 1] T . The deterministic equations can be 
derived from this master equation and yield: 

^ = J2rW(X dcU r) 

r 

The Wentzell-Friedlin theory [5J [B] assumes that for any path [X] in configuration space: 

P([X])~cxp(-S([X])) 

with <S([X*]) = for the deterministic path [X*]. It follows that the dominant contribution to the 
probability of a fluctuation that brings the system to state X starting from a stable state X m can be 
written as: 

P(X|X m ,t = 0)=exp(-S(X)) (7) 

where 

S(X) = min [X ] : x m -»x<S([x]) (8) 
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where the minimization is over all paths [X] starting at X m and ending at X. For X far away from the 
steady state, the probability of occupation P(X) is equivalent to logarithmic accuracy to the probability 
of the most likely fluctuation, V(X\X m , t = 0) that brings the system to X. The assumption of the form 
given by Eq. [7] for the occupation probability is known as the eikonal approximation. 

Using a smoothness assumption for W(X, r), and since the changes in numbers of A and B nodes 
are O(l), we can neglect the difference between W{X — r, r) and W(X,r). With this approximation, 
the eikonal form for the occupation probabilities in the master equation yields the following equation for 
SCX) i: 

ds 



^x,-j=0 (9) 
where 

(x, p) = ^(*, r)(cxp(rp) - 1) (10) 

r 

and 

x = X/N, w(xl, r) = W(X, r)/N, s(x) = S(X)/N. 

Eq. [5J is analogous to a Hamilton- Jacobi equation for the action of a system with Hamiltonian given by 
Eq. 1101 The corresponding Hamilton equations of motion for components of position x and momentum 
p are: 

Xi = ^— Pi = (11) 

dpi Xi 

with s(x) playing the role of the classical action: 



([x]) = / L(x, p)dx = / pxdx 



'[X 

where [x] denotes a particular path obeying the equations of motion fEqs. [TTj). 

Following this Hamiltonian formulation to characterize the fluctuational paths of the system, our goal 
is to find the path with minimum action that reaches the separatrix in phase space of the deterministic 
motion, starting from the vicinity of the stable state under consideration [6l [4]. Arguments in [3] show 
that the fluctuational path reaching the separatrix with the minimal value of the action, is the path 
that passes through the saddle point. This is the optimal escape path, i.e., the path whose probability of 
occurrence dominates the probability of escape and we denote it by [x opt ]. This path can be found by 
integrating the equations of motion Eq. II H and finding the required path that starts from the vicinity 
x m to the saddle point x sac idic- Thus following Eqs. [3 [8] we have for the probability of escape from the 
current stable point in which the system is trapped: 

-Pcscapc = -P(x S addlc) ~ CXp[-iVs(x sadd i c )] (12) 

where 

s(x saddlc ) = / L(x, p)dx 

•'[Xopt] 

and the transition time (or time to escape from the steady state) follows: 

Switching ~ exp[iVs(x sadd i c )] (13) 

In practice we start from some point x in the vicinity of the stable state, and to obtain the corre- 
sponding momenta p and action s(x), we employ a Gaussian approximation [4]: 

S{*)=Y J Z lJ {x l -x?)(x 3 -x™) 
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where Z satisfies an algebraic Ricatti equation: 

QZ- 1 + Z- 1 Q T + K = 

where Q is the linear stability matrix (Eq. [3} evaluated at x TO and Kij = J^ r , u)(x m ,r)r,-r J -. Solving this 
Ricatti equation yields Z which in turn yields S'(x) and p(x). 

In order to find the optimal fluctuational path of escape from a given steady state, we numerically 
generate fluctuational paths from various points close to the steady state (we explore points at intervals 
of 10 -5 along the x\ dimension and 10~ 2 along the X2 dimension around the steady state) and find 
one that passes close enough (no greater than a distance of 10~ 5 ) to the saddle point. The equations 
of motion, Eqs. 111! are integrated using a trapezoidal rule to generate these paths starting with initial 
conditions obtained using the Gaussian approximation described above and subsequent numerical solution 
of the Ricatti equation (we use a Matlab Ricatti equation solver for the latter). The scaling behavior 
of switching times obtained using this approach for various committed fraction values as a function of 
distance from second-order transition (or cusp) point are shown in Fig. 5 of the main text. 
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